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I'NITED  KIXCDOM 

ABSTRACT 

The  instability  of  a  viscous  fluid  inside  a  rectangular  tank  oscillating  abiuii  asi  axis 
parallel  to  the  largest  face  of  the  tank  is  investigated  in  the  lin<'ar  rx'gitne.  TIk'  Ikju  is  sIkmvu 
to  be  unstable  to  both  longitudinal  roll  and  standing  wave  instabilil i('s.  The  particular 
cases  of  low  and  high  oscillation  frequencies  are  discussed  in  detail  ami  the  ix’sults  obtained 
for  the  standing  wave  instability  at  low  frequencies  shed  light  on  tlu>  corrt'sijondiiig  steady 
flow  instability  problem.  The  relationship  between  tlx*  roll  instability  and  con\'('ctive  or 
centrifugal  instabilities  in  unsteady  boundary  layers  is  discussed.  TIu'  ('igruivalue  ))roblems 
associated  with  the  roll  and  standing  wave  instabilities  are  solved  using  I'locjuet  theory  and 
a  combination  of  numerical  and  asymptotic  methods.  The  results  obt aim'd  are  comi)ared 
to  the  recent  experimental  investigation  of  F3olton  and  Maun’rf  1992)  v\liiui  indeed  |)ioviiied 
the  stimulus  for  the  present  investigation. 
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1  Introduction 


Our  cunct'ni  is  wit  h  ilu'  iiist ahility  i)f  I  lie  tiiistt-aiiy  lluw  of  ,1  vi-cous  Ihijd  in-ido  .1 
rcct angular  tank  oscil !at  mg  al)out  an  axis  paialinl  t o  oiir  fare  (d  i  Im'  t  an io  I  ho  orient  ;i!  ion 
of  t!i('  axis  ol  oscillation  is  shown  in  Figure  (I.!).  An  ex[H'i  inieiit  al  iineOigaiiou  o!  tlie 
diffeiauit  reginit's  tor  tlu'  flow  insidf'  tlie  tank  is  <les<'ril>ed  in  iFilton  ,ind  Maurer  t  I't'd'Jt. 
heralter  referrt'd  to  as  [-{.M.  If  t  lu'  g<'om<'try  ol  the  tank  shown  in  Figure  1  i.l  i  i--  -m  ii 
that 

-i  ->>  1  —  i 

,/  '  ,1 

tluMi  the  basic  flow  is  unidirectional  (in  the  .c"  dir<-ctionl  ami  ile[jendent  oiib.  on  //" 
atid  lime.  In  tact  tlie  flow  is  an  <'xact  solution  of  the  .Xavier  Stokes  e<|uaiio!,  .  Ih.e 
stability  of  the  flow  is  therefore  govenuHl  Iw  a  system  ol  partial  dilferentjai  erin.itious 
witli  coefficients  periodic  in  time.  We  shall  show  tiiat  the  linear  iusi,ibility  j)rob!<’m  i' 
closely  related  to  tho.se  for  r('ntrifugal  and  coiivectiv<>  insi  abilit  ies  in  I  iine  perioiiic  flows. 
Before  discussing  the  results  cjf  H.M  we  shall  first  briefly  re\i(nv  the  iclewtnt  details  of 
the  related  flows  mentioned  above. 

Fhe  experimental  and  theoretical  invest igaticuis  of  tlie  centrifugal  instability  of  a 
Stokes  layer  by  Seminava  and  Hall  (  l9T(i)  showed  that  time-periodic  llovvs  can  sustain 
instabilities  not  acressilile  to  quasi-steady  instaliility  I lit'orii's  based  on  the  insi ant aneons 
velocity  profiles.  Seminara  and  Hall  (H)7b)  found  tliat  a  torsionally  osi  illating  crdimlct 
in  a  viscous  fluid  drives  an  unsteady  houndary  layer  which  is  mislaiile  to  davliu'  vortex 
like  instabilities  at  high  enough  frequencies  of  oseillation.  .A  more  del  ailed  expiu  imeiital 
investigation  of  the  problem  by  Park,  Barenghi  and  Donnelly  (1!)S())  conlirnu'd  the  s('c- 
ondary  .subharinonic  destabilization  of  tlu*  most  dangerous  mo<l<‘  found  by  Seminara  and 
Hall  ( 1976).  An  approximate  descriptioti  of  this  subharmonic  Iweakdown  wa--  later  .given 
by  Hall  (1981).  Subsequently  it  was  shown  hy  Hall  (l‘)81)  and  Papageorgin  (PtS7)  that 
the  instability  mechanism  found  by  Seminara  and  Hall  can  occur  in  spati-dly  localizetl 
positions  in  more  complicated  unsteady  boundary  layi'r  flows. 

The  convection  mode  of  instability  of  tlie  unsl(>a<iy  thermal  houndary  lavcr  in  a  sc'uij 
infinite  mass  of  fluid  adjacent  to  a  tiim  periodically  la’aliai  rigid  wall  was  invest  igati-d 
theoretically  by  Hall  (198u).  In  tlie  latter  jiaper  it  was  shown  that  tlie  cunrcci  loii 
problem  for  a  fluid  with  Prandtl  numlier  eipial  to  unity  is  idiuitica!  to  that  gova-rniiig 
the  centrifugal  instability  [uo!  Icn.  6/.  t.i.  fluv.  io  n  oiphus  rotating  cvlmdcr 

in  a  uniform  stream.  As  y<’t  only  the  .severely  trum  aletl  e(|uatioiis  disi  usscd  b\  Hal! 
have  lieen  solvx'd  and  the  results  found  suggest  that  lli<-  most  dangerous  mode  a 


su!)ltar!iu)nic  uiu'.  If  tliat  is  llu-  only  uiistahU'  iiiudf  llunt  iinll  >  analysis  siiy<4i-~ts  tha! 
the  flow  oil  a  ra[)i<liy  rotaliiisj;  cyliiHliT  is  staMc  siiiro  in  ihai  jin^iili'in  tin-  [xflar  anyh' 
0  [)lays  t.h('  I'olc  of  tiuK'  and  solutions  jH-riodic  in  0  with  [ift'io*!  an-  ot  no  jdiV'iiai 
n'levaiu't'. 

Allot  lior  con  vt'i't  ion  [iroldcin  assixialrd  with  1  iiiif-|M‘iiodic  lorciny  i-'  that  disiiis'fd 
hy  (ii'cslio  and  Sani  (l!)7ll).  I  hc  latter  antliors  investigated  the  instalidity  ol  a  layc’ ol 
fluid  heated  stt'adily  IVorn  lielow  in  a  t  iine-pc'iiodic  gravity  held.  It  wa,>  lound  that  the 
staliililv  probh'in  is  govt'ined  lyv  Mathitni's  etpial ion  atni  t  hat  the  donnnant  lns1al)iht\  is 
snhharinonir  one.  hhe  |)rol)leiti  discussed  by  (Iresho  and  Sani  is  (j|  considerable  pract  ical 
iin|)ortance  becanst'  of  the  presence  of  convection  in  a  inicro-yravity  ('iiviroinent  when- 
viltrations  cause  tlie  (dfet  tive  gravitat ional  liehl  to  be  oscillatory  in  time. 

I  he  possible'  instability  of  tini<'-i)e'riodi<'  Hows  to  I  rax'elline;  wa\'e  disi  nrbaiK  <-s  has 
l)y  contrast  not  rece'ived  much  atte'ntion.  The  fnndatm'iital  pretldein  In're  concerns  t he 
liiK'ar  instability  of  a  Stokes  layt'r  ott  a  traiisve'tse'ly  oscillating  ri.yid  plane  w;dl  to  waves 
jn'opagating  in  the  flow  direction,  dliis  prolrh'in  was  first  investigated  by  Kerezek  and 
Davis  (107!)  who  found  that,  even  though  the  instantaiie'ons  |)rolih’s  can  Ik'  highly  in- 
fiectional  and  therefore  massively  unstalile.  the  oscillatory  flow  between  parallel  [tlales 
is  stable  according  to  a  Flocpiet  approach.  Later  Mall  (107iS)  showt'el  tliat.  even  though 
the  Floquet  solutions  of  Kerezek  and  Davis  were  greatly  d('p('nd('nt  (ui  tlu'  presence' 
of  a  stationary  wall,  the  Stokes  layer  ott  a  wall  oscillating  in  a  viscous  fluid  is  stable. 
However  tlie  instantaneons  velocity  profiles  associated  with  a  Stokes  layer  can  Ik'  mas¬ 
sively  unstable  on  the  basis  of  a  <inasi-steaily  analysis.  Such  an  approach  is  \alid  at 
large  Reynolds  numbers,  howe'V'er  the  unstable'  .solutions  canne)t  lie  e  enitinue'd  ove'r  a  full 
perioel  of  ejscillatiejM  to  produce  Floepiet  solutions.  More'  re'cent  work  by  .Akhavan  e't  al 
{ 1991a. b).  based  on  full  numerical  simulations  of  the'  Navie-r-Stoke's  e-epial iems.  sugge'sts 
that  transition  to  turbide'uce  in  Stokc's  layers  can  be  attributed  to  higher  ejrde-r  instaliii- 
itie's  asseiciateel  with  the  jmimary  instaViilitie's  of  the  inst antane'emsly  invise  idlv  unstable' 
ve’leK'ity  pre)file  of  the  basic  state'. 

We'  shall  nenv  discuss  the'  main  re'sults  found  by  HM  in  the'ir  e'xpe'iime'nt al  inve'sliga- 
tion  e)f  the  flow  in  a  flapping  rc'ctangular  tank.  A  me>re  eh'taih'd  eliscnssioii  e>f  the'  re'sults 
can  be'  feim'd  in  (;4  of  this  pa|K'r. 

In  oreler  tei  characteri/.e'  the  freriiioi,ry  ‘»f  tie'  fhc.v  DM  nP ''e>duc('<i  tin'  parame-t-'- 
(p  ~  wlie'ie'  tv  is  the'  fre'e|ue'ne’y  of  ose  illation,  the'  kine'inatic  viseosity  ami  il  as 
shown  in  Figure*  (1.1).  At  a  fixe-d  value'  of  the-  flapping  angle*  o  HM  eibse'rved  that  for 

•) 


small  t'uoiigh  \aliu's  ol  tlu*  fl(j\v  was  slahU-.  W  lu*ii  4>  w;is  inriasiscfl  a  liitiu  <  at  ioii  tu  a 
weak  roll  stat<’  took  pkua'  at  a  critical  \’aliie  ot  <l>.  II  this  critical  \aliie  el  ‘k  is  (h'lioted 
l)y  then  BM  sliowed  that  <I>,  j  is  a  moiiotonically  decreasiiig  functioii  iif  o,  At  a 
second  critical  valnt'  of  a.  a  strong  roll  state  was  lound  h\'  BM  and  ;.t  suflicient  l\ 
smalt  values  of  a  this  mode  exhihited  hysteresis.  At  higher  values  of  d>  was  \  modes  v'o  re 
observed  experimentally  though  BM  suggest  these  w<'re  associated  with  (Uid  elfects.  At 
very  high  values  of  a  turbideut  How  superiinpo.sed  on  some  la-sidiial  rt,!e  stun  lure  was 
observed  in  the  e.xperiments,  I  h<'  analysis  in  this  j)ap<’r  wdl  hjcus  ois  the  origin  of  the 
strong  vortex  state  found  experimentally.  how('v<'r  our  results  will  also  suggest  a  likeb 
candidate  for  a  mode  responsible  for  the  onset  of  the  wavy  states. 

In  the  following  section  we  shall  tormulate  the  liin'ar  insiabilit\'  probietn  lor  tlie 
unidirectional  How  in  an  infinitely  long  Happing  tank.  The  ('fptat  ions  we  diui  ve  govi'rn  the 
linear  instability  of  the  How  to  disturhanct'  periodic  in  th(>  .r'  and  c"  direct  ions.  In  Sect  ion 
3  these  equations  are  discussed  for  roll  modes  which  are  laktni  to  he  indojiendent  of  .r\ 
The  particular  cases  of  large  and  small  4>  are  disenssed  in  ^3  whilst  in  1  nuiiK'rical  o'sults 
for  4*  of  0(1)  size  are  presented  atid  onr  results  eompaixul  to  (wix'iiment al  observations. 
In  Section  5  we  disruss  travelling  wave  clisturltances  which  ar(>  iinh'itetident  of h’inally 
ill  §')  we  draw  some  conclusions. 

2  Basic  flow  and  the  stability  problem 

Consider  the  flow  of  a  viscous  fluid  iii  the  rectangular  container  deliiied  hv 

-  <  X  <  --  <  11  <  -A,  <  <  A,  (d.l ) 

with  respect  to  a  Cartesian  coordinate  system  (.r*.  i/',  j:*  ).  The  fluid  is  takt-n  to  he 
incompressible  and  the  density  aiifl  viscosity  are  ilenoted  by  p  and  n  res|)('ct  i  vel\'.  I  he 
fluid  is  set  in  motion  by  the  oscillation  of  the  container  about  the  axis  with  angular 
velocity  (0,  0,  sin  cnC).  Following  B.M  we  define  the  fr('(|iiency  parameter  4>  by 

4>  =  - — 

i' 

so  that  4)  >  1  corresponds  to  a  situation  where  viscous  ('ffects  atx'  small  wliilst  4>  <  1 
corresponds  to  a  viscous  dominated  flow.  W<‘  can  suppose  that  I  lie  vtdocit  y  and  pressuK- 
of  tlie  fluid  are  scaled  on  aud  ami  opu>^d^  res|)ectiv('ly  whilst  dimensionless  varialiles 
ix,y,z)  and  I  are  defined  by 

(.r,y,^)  =  /,--).  t  ^  ^l\  (A3) 
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vvlit'i'f'  /*  tlt'iK)t('s  tiuK'.  W’itli  rcspoft  to  t lio  ruordiiiatr  system  iiKA  iii!.’  wi'li  !!!<  tank  the 
Navii  r  Stt)k('s  e(|uatiui)s  take'  the  form 


(livxi  —  0. 


u,  +o(u  ■  V)u  — -All  =  —  V/id- 2o  sin  / 
<n 


wliirh  must  be  solved  subji'ct  to 


■A 

( 

—  ti  j 

+  a  sin“  f 

-f"  ( ’os  / 

“  ./■ 

V  ) 

\  '» / 

(2,1) 


u  =  0.  cju  j'  =  i/.pd  *.  jj  ~  i  — .  c  —  ztl.-il 


(2.0) 


If  we  write 


•siif^/,  >  ■>, 

p  =  -  {.r  +  //  )  -  T)j  eos  t  4-  p 


then  (2.1)  IjeroiiK'S 


( 

A 

^  2y  cos  /  '' 

ti,  +  o(u  •  V)u  —  =  —Vf)  +  2a  sin  t 

-u 

+ 

0 

V  0  i 

and  for  eonvenienee  we  now  drop  the  A  notation. 

In  order  to  make  analjdica!  progress  we  assume  that  Lj^/d.  L:/d  are  large'  so  that  we 
ran  firop  the  lioundary  conditions  at  x  =  dzLrd~K  z  =  ±L;d~^.  As  in  BM  this  enables 
us  to  look  for  a  unidirectional  flow  of  the  form 


u  -=  (iiiy.  f ),  0, 0), 


V  —  w  =  0,  p  ~  —2a  sin  f.  /  udy. 


vvliere' 

—  +  Ht  =  2y  cos  f . 

a  =  0.  y  =  ±7y 

I’he  r('<iuired  solution  is 


(2.7) 


wlu're 


«  =  +  Q}(‘'  +  COMPLKX  (’0N.B;C;A'I'K, 


Q  = 


i  siidi 


2sinh 


12.S) 

(2.<)) 
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I'or  large  valiu'.s  of  0  tlu'  iini(  tio:i  ()  is  ex[)uiienl  lalty  small  away  lioiii  layers  oj  deptli 
0(0 )“ 5  iK'iU'  i)  —  iy.  rims  for  large  0  we  liiid  that  nc-ar  //  =  ■; 


I  ''0  1  I  '■  I 

II  =  sin  /  —  sin  /  +  \/  — ( 7 - )  m  -  ■ 

\  2  '  2  ' 


Kor  small  vahu's  of  0  ihe  llnid  vespoinls  iji  a  <inasi-s1('a<ly  manner  1o  lli<‘  Itaclng  am!  we 


ol)tain 


cos  '  m  .  .  ■ 


i'O.l  !  I 


riie  flow  in  a  tank  rotating  with  constant  angular  xehxity  is  then  obtained  hom  the 
al)ove  i)v  setting  /  e(|ual  to  zero. 

We  now  perturb  the  basic  flow  given  by  (2.iS)  by  writing 

u  =  (u.0.0)  +  [U(/y.  /)exp?  {A.r  +  k-}  +  ('OMPLhLX  ( ’().\.) Ct bVl'K]  +  .  .  . .  (2.12) 

If  w('  assume  that  jUj  <C  |u|  then  vve  can  linearize  tlu'  ecpiations  (2.(i)  tc^  give 

l\('  +  V„  +  ikW  ^  0. 


Cl-  +  a.Vii,i 

cv 

CM- 

Here  the  operator  C  is  d<'fitied  by 


=  —i\P  +  2(\Vs\nf. 
=  — /y  —  2f>/ '  sin /. 

=  -ihP. 


C  =  -[M  -  A"  -  P]0-'  +  InXU  +  (),. 


(2.11) 


rhe  eriuations  (2.12)  must  Ire  solved  subject  to 


n  =  v  =  w  =  o,  !j  =  ±- 


(2.I.V) 


Since  «  is  a 


be  found  with 


function  of  time  we  anticipate  that  solutions  of  (2.12),  (2.1”))  may 


{(I,  v;  VV',P)  -  r'‘'(fMA  W.  P) 


(2.1(.) 


where  lf.V,W  and  P  are  periodic  with  respect  ter  /  and  tin'  I'lofpH't  e.xpoiK'iit  //  is 
complex  and  is  a  function  of  0,  rr.  A,  k.  d'lie  stability  of  the  flow  is  then  determined 
by  the  sign  of  //-r;  if  solutions  of  th<'  form  (2.  lb)  e.xist  with  //,  >  (1  then  the  flow  is 
unstable.  In  the  next  section  we  discu.ss  solutions  of  the  eigenvabu'  piarbh-m  ior  the  <  ase 
A  =  0  which,  following  H.M  ,  we  refer  trr  as  “roir'  modes.  In  se<  tion  ■”>  wc  inv('sliga!c  the 
possibility  of  Tollmien-Schlichting  wava-  instalrilitir's. 


3  Roll  modes  of  instability 


We  shall  now  seek  solutions  of  (2.15)  with  A  =  0.  it  is  llu'U  ronvenien!  to  eliminate 

ir  am!  P  to  give  the  iollowing  eoupled  pair  of  equations  for  /  ami  1: 

-  ctppp  =  «ci>v  . 

{0l  -  -  k^W  =  -2a<^k-s\u  I.  (.Il  l 

P==V'  =  V,  =0.  = 

It  is  oI  interest  to  note  that  (5.1)  al.so  govern.s  the  stability  of  a  vertically  oscilhit  iiui, 
Boussinesq  fluid  between  parallel  walls  y  =  ±^.  hi  that  case  the  Ibiid  has  Brandt! 
number  unity  and  the  upper  and  lower  walls  have  temperature  pro|)ortional  to  isiii/ 
respectively  whilst  the  fluid  is  subject  to  a  gravitational  field  proport iotial  to  sin  t.  In 
that  case  the  Rayleigh  number  for  the  flow  is  so  that  at  0(  1 )  valties  of  <l>  we  should 
anticipate  tmstable  solutions  of  (3.1)  for  o  =  0(1).  Before  discussing  the  numerical 
solution  of  (3.1)  for  0(1)  values  of  a  and  <I>  it  is  instructive  for  us  to  first  consider  thi’ 
further  limits  <f>  -e  0  atid  4>  oc. 


Low  frequency  limit 


Iti  the  low  frequency  limit  (3.1)  redttces  to 


(Of  —  k^  —  ^0t}(’  =  — 2('v<I)  <  silt  /  +  <f>  cos  t 


r 

0 


_1_ 

2A 


<1>^  sin  t 

12 


{c)l-k^ 


k^)V  =  -2rvtI)/('V'sin/.. 


2  210^  j 

(3.2a) 

i3.2b) 


f '  =  P  =  W 


U  =  ±;^ 


(3.2c) 


WV  now  indicate  how  a  W  K  H  type  of  solution  of  the  aliove  equal  iotis  can  Ih'  found.  If 
we  set  rt  =  0  tlien  tht'  ('quations  for  V  and  V'  decouple  and  it  is  easy  to  see  that  the  flow 
is  stable  with  decay  rates  iif  size  on  the  /  timescale.  W'e  anticipale  that  this  ih'cay 
will  theti  be  balanced  by  grovvtli  associated  with  the  ajiparently  desfabilizitig  t.ertns  on 
the  right  hand  sid<'  of  (3.2a.b).  rids  is  achiev(>d  if  o<I>  ~  0(1)  so  tiiat  we  write 


and  lei  (/".  V')  =  o(/y)- ^()(//)) ^  j' ^*1^}  +  ••••  I  Im  eiiyenvaluc  prulilrm  tni  ijic 

local  growth  rate  rr  is  then  found  to  lie  givtui  !)y 

l<il  —  A-"'  -  (t]/ u  =  -2.  l!?iii  /\  o. 

A  i)  =  h()  =  1(1^  !/  =■  i-. 

Thus  t  appears  otdy  as  a  parameter  in  the  z<'roth  ord<’r  i)rol)l<‘m  so  we  have  an  oi(litiai  \ 
differential  system  to  determine  the  eigenvalues  rr  =  a{  I .  J\ .  A).  In  fa<  t  hy  ii  plarinu 
V’osint  hy  Vq  we  see  that  (3.d)  is  then  ecpiivalent  to  the  Henard  j)rohlem  hu  a  iluid  ot 
Prandtl  number  unity  between  the  plate  n  =  ~\  at  tempcuaiure  units  and  the  plaif 
t/  =  ^  at  zero  temperature.  However  the  effoctiv<>  Rayleigh  nmnlK'r  is  —  1.1*  sin"  /  so  t  hat 
the  flow  is  stable  and  rr  must  have  negative  real  part.  This  I’an  b(’  seen  trom  f'b.d)  for 
large  values  of  .4  by  writing 

<7  =  (TuA  + - 

f  0  =  (  00  +  .  •  .  -  id.  I  I 

ho  =  h'oo  +  .  . . . 

The  eigenvalue  problem  for  rrg  then  becomes 


'7oK-A-^]Vo  =  U-'sin'/Vo 


Id-'.) 


h'o  =  0,  y  =  ±:^. 

We  notice  here  that  the  limit  .4  — >  oc  is  an  inviscid  one  so  that  the  eigenvalue  problem 
is  now  associated  with  a  second  order  differential  equation.  In  fart  tlu'  eigcuivahu's  of 
(3.5)  are 

±2A:r  sin  t 

V,.v  +  e’  =  ''‘"I 

which  means  that  the  flow  is  inviscidly  stable.  Thus  we  have'  shown  abov('  that  for  4>  ;>>  1 
a4>  ;§>  1  the  disturbance  has  an  imaginary  growth  rate  of  size  0(o). 

The  above  discussion  shows  that  no  neutral  disturbances  exist  for  r><h  =  ()( 1 ).  T  -C  1 . 
In  order  to  obtain  neutral  disturbances  we  must  increase  o  until  suflicient  ('xponent  i.d 
growth  takes  place  near  the  times  wlien  sin/  =  0  to  l)alance  th('  ('Xjxuieiitial  decay 
associated  with  viscous  effects  at  other  times,  'I'he  exponential  growth  takes  phu o  in 
0(<hj  time  intervals  near  the  times  when  sin  /  =  0;  this  nu'ans  that  a  disturl.ance  grows 
by  an  amount  of  order  for  some  constant  (A  Viscous  efh'cts  on  the  other  hand 


Icail  to  (lcca\  !)>■  la(fors  size  r  ^  ^  for  somo  f),  I  Inis  the  raft'  (!<•<  m.ist’s  wifh 

ii  o  atitl  4'  art'  held  fixt’tl.  llowt-vcr  when  k‘  unrnascs  to  l)(  4>  I'  ■  \iscou>  “flfits  in  iIk* 
l.ulk  (.1  [  tilt'  Ilow  al  so  foiiit’  into  play  ainl  a  iiiiiiinHiin  rate-  is  a<  Wo  1  licrolttif  si  rk 

a  Solution  ol  id.L’)  ior  the  small  4>  cast'  with 

Ii  “  ^(<{>  ' .  !  d. 7<i ! 

k-  =  A  'l>- ■ 

lit'itl  iixod. 

\\('  then  ('Xjiaml  /  ami  I  in  the  form 

{l\  \")  =  }.  \„{n.l))<\'^vxv 

r,=^) 

v.  il  h  ( ’i  =  I =  (Ti  =  0. 

i  lit'  leading  order  systems  lo  tlett'niiim' tT] .  tx.)  art'  then  ionnd  to  Ix' 

tx,,/ D  =  2s'\ul\uli 

mi^i)  =  -'isintf  id.'.)) 

(•■{.10) 

{••{.in 


(d.rj) 

(.{.!:{ I 


(X„ /  J  —  'ls\ut\  [>fi  —  fT^i  'll 

iT,|l  j  =  -'i.sill//  2^^  ~  ^2^0  + 


ttrii  1 .1,,,^ 
7- 


1  he  ( (insistency  of  (■{.0)  retpiirt's 


tx,";  =  —  1  siii^  /  li- 


so  that  tXo  is  purely  imaginary  and  th<>n  (-'{.lO)  is  consistt'nt  of 


{  llvv  ■”  —  0 

^11 


and  the  solution  ot  this  ('(|uat,ioii  which  vanishes  at  ±7  ~  /;  is 


It)  —  sin  nrrfi/  -p  —  ), 


With  f'(l  )  to  1)('  determined  and 


nil< 


Vi 


It  --  i .  2.  d, . . . 


4>-’d  1 

I  YirrjD^iji 

J 

t,;=U 

s 


i  has  rr-j  is  alsi;  pun'ly  imaginary  and  iiidard  (Tr,.  nr.,,  rrs  arc  alsn  imaginai  v.  1  he  c<|i;al  ii iii'' 
to  (Icterminc  <t,;  art'  found  to  he 


(T()f  V;  -  ~B  sin/ Vi 


fToVis  -t-  2B  siri/r,; 


-  [^2'  1  d-  ,(  +  1  +  ('^t;  +  )l  n 

+  'lB\(j\  eos/l  {  +  '/2  •■'ill/ 1 ’2] 

—  [<72  1  il  +  (T:\  Va  +  t7.,  I  2  +  tT;;  I  i  +  (  7,,  -f  K']\  f,  j 

+  A'~‘[7,i\it  <^A  1  -t-  7;, I  1  +  7.ili, 


whe 


re 


<l\  = 


<ii  = 


^-4--dL\/|2 

2  2 !()  /  ’ 


CM  1 


CMC 


2  2d 

it  we  t'liminate  /  0.  Vq  ironi  (0  ,•.}  we  obtain  a  differt'iit  ial  t'fjnatitni  for  V,.  An  examina 
tioi)  ot  the  (list nrlrance  structure  in  the  viscmis  wall  layers  shows  that  tlx*  etpiation  lor 
V,|  must  Ire  solved  subject  to 


I  >  n<-(-  1 )" 

=  ~r  •  n 


w 


CMb) 


//  =  -7- 


If  the  e<|uati(rii  for  V't  is  to  havt'  a  solution  satisfying  thest'  condititms  tlit'u  t  lit'  real  [rart 
of  <7,i  is  given  Iry 

I  sitd|i 


l<^ 


-  K\ 


(2.17) 


riie  amplitude  function  C{i)  is  determined  at  higher  order  and  is  found  to  bt'  singular  at 
the  instants  wlien  7o  vanishes.  In  order  to  find  the  disturlranct'  structur<'  at  su<'h  tinu's 
we  consider  a  small  time  interval  near  for  example  /  =  0  and  deiint' 


T  =  4)-'/. 


(:Ms) 


Phis  scaling  is  implied  by  {3.2aj  which  sliows  that  the  first  and  second  ti'it'is  in  th<' 
brackets  on  the  right  hand  side  of  this  erination  are  comparable  whemwer  sin/  ~ 

We  then  note  that  (2.2)  may  then  be  written  in  the  form 


~~dr)B 

4)T 


■2B 


I  \2 


4) 


y  .  i  +  Mii  tiii! 

12  2  '  2 


<t>  ^ 


'y  -  V 


2Bh 


t.'Mlta! 


id.l'tS.t 


\  'P '  /  \  ‘P  <  /  ■■ 

1  hese  erpiaticrns  may  be  soivc'd  usitig  a  WKIf  api)roa<  h  to  take  l  ari'  of  the  time  de|)en 
deuce  and  by  noting  that  for  small  valut's  of  (y  —  “)  wr  can  i'2,nore  the  t  Idl'd  and  foiirlh 


tcnns  ill  the  lirackt't  nii  tlic  liainl  side  ul  (ii.nta).  In  oralci'  to  halam c  tin*  lime 

dc'rtvat  iv('  with  tonus  on  t  ho  i'itj,li1  haiul  sido  ot  (-i.!’))  wo  must  tako  ()/  -v  uni) i  and 

tiio  fi  (lono'uloiH’o  of  I  ho  (list tirl>a!Ha’  thou  shiiiiks  to  a  thin  la\'or  (h  ihi(‘knt's>  <!>■■  nc.ii 
tj  =  i.  Ill  hid  a  similar  layer  exists  near  i)  ~  —  ~  Iml  the  slrudiiro  is  similar  fo  that  at 

the  njipor  wall.  Wo  (h'fiiio  “  (//  --  and  ihoii  look  for  a  solution  oi  (d.l'.ii  luaii 

//  =  X  ^»t  t  lit'  lorm 

ilW)=  id-ifiii 


■■ 

(/  n(^.7  ).lu(^.7  )  4-  (f)A/  i(^,7').  l,(^7'))...j  exp 

/  iUT) 

4-  <l>A/,l7  1  -  . 

.  ,  uil 

(p  -  ./ 

, 

If  W('  snhstiluto  tlu'  above  o.xpaiisions  into  (.'hi!*)  and  O(|nato  lorin.s  ol  the  [a.'Wi-rs  in  <I> 
wo  obtain  at  xorotli  order  a  pair  of  linear  otiualions  lor  I  \(i.  1  bo  consist oncv  ot  those 
('fpiat ions  yields 

.y  =  -i/;-r(r+^)  loai 

so  that  wo  have  an  oNjnmontially  ^vowiii"  stilution  in 


•—  <  T  <  0. 

!•> 


\V(’  a.ssiimo  that  I  is  in  this  ran»,o  and  o<j,isidor  the  root  of  (d.2i  )  with  ./o  >  0.  .At  no.xi 
order  w<'  (iiid  that  the  linear  ('(|uation  for  l\.  ly  obtained  are  eonsislent  il 


<p\;>  2,/,,..,.  ,  i-iA'-’^,, 


<ie 


(d.22) 


I'his  ('(piation  is  then  soh'od  subjerf  to  Vh  =  0.  <^  =  0  and  such  that  l„  —»  0. 
riiis  enabh's  us  to  express  Ij,  in  terms  of  solnlions  of  Airy's  ('(piation  and  the  (juantity 
./i  can  then  la'  exjiressed  in  terms  of  the  zero's  of  .1,.  1  he  solution  (d.22)  lails  when 

7’  =  0.  7  =  — p  and  WKB  turning  point  layers  (with  la'spe'rt  to  /)  art'  tu'eded  to 
contH’ct  (.’hS)  and  (.■}.22).  .Across  tlic'se  layi'rs  Uie  two  oscillating  solutions  (2.8)  witli 
<T(,  =  ±2(sin//7  conned  with  the  ('Xponi'iit iaily  di’caving  and  growing  solutions  ('hdO) 
with  •/((  =  ±2/:?[  — 7  (7  +  A  (leriodic  .solution  is  obtaiiu'd  by  choosing  U  such  that 

t  h('  ('xpoiK'iit  ia!  growth  in  —  dj  <  7’  <  0  is  identical  to  tlu*  ('Xjionent  iai  (h'cay  associati'd 
witli  (T,;  in  (2.8).  VV('  not('  lu'rt'  that  tlu*  |)articular  form  ol  tin'  tinu'  d('|)('nd('nc('  ol  (2.!*) 
enables  us  to  consider  only  tlie  int«'rval  1)  <  I  <  it.  If  wi'  tln'ii  considf'r  the  h’ast  decay  ing 
solution  (2.17)  with  n  ~  !  we  find  that  the  smallest  value  of  li  which  leads  to  a  neutral 
solution  of  (2.d)  lor  I  satisfies 


XX  /■  |sin/!^(// +  A'‘;r  ^  2/y  r  [ -7'(7’ +  )]f 

l  \  ^  Jn  7-  A  12 


(2.22) 


10 


"Fhe  al)Ove  ('(juatioii  can  then  he  seiveii  lor  li-  (iinplicilly  asMinied  to  In-  po.jm-  in  iln’ 
tlerivation  of  (.‘hlh'}))  a.s  a  fniiction  of  A’.  Fioiin-  (■i.I)  shows  H  as  a  ol  A.  '.ve 

see  that  B  ~  A,"^  B  ~  A,~‘  for  small  and  larjte  A  resperi  i  vely  and  that  B  aiiaiiis  ,i 
minimnm  at  some  intermediate  value  of  A  . 

If  the  integral  on  the  left  hand  si<le  oi  (.'{.2d)  is  integrated  nnmeiically  we  find  that  the 
minimum  occurs  when 

B  =  7Mm).  id.dlai 

A'  =  7.!)9.  (.{.'JUj! 


Thus  tlie  most  dangerous  mo<ie  for  1  has  o  given  liv 

Tdhdo 


ft 


<p  i 


+ . . . 


We  postpone  fttrtlier  discussion  of  (i{.'2u)  until  the  next  section  wluu'e  we  discuss  i  ht- 
numerical  solution  of  the  eigen vahu'  problem  for  4>  =  (J(  1 ). 


The  high  frequency  limit 

For  large  values  of  the  freciuency  parameter  the  fum  tion  ((^(// )  appearing  in  t  he  si  ability 
equations  (3.1)  develops  boundary  layers  of  thickness  ^"5  at  //  —  and  is  <'xpon(':!t  ially 
small  elsewhere.  It  follows  that  any  instability  must  be  locali^ted  in  tlu'se  layers  Pjr 
definiteness  we  focus  on  the  layer  at  t/  =  and  define 


(3.'2f.) 


The  dominant  terms  on  the  right  ami  left  hand  sides  of  (3.1)  then  balance  (or  ?/  =  ()( i  j 
if 

(  ~  O'  \  4*  2  _  \  ~  o  f 

with  k  ~  O)*!*?).  Hence  we  must  take  o  ~  ami  write 

B 


o 


+  . . . 


+  ... 


k 


aiu!  the  zt'i'oth  ordt'r  approxiinatimi  tu  (d.l)  in  llin  lower  wall  layer  can  he  written  in 
the  fonn 


(t);;  -  4-  COMfM.KX  ('()N.ir(iATK| 

(t^;;  -  -  Ot)i()f,  -  A-‘)K  =  -2Bk^s\utr. 

/  ■  —  V’  =  \  =  0.  >)  —  0.  /  ’.  V'  0.  j/  — +  oc. 

Solutions  of  this  system  of  the  form 


ran  be  found  and  the  Floquet  exponent  //  is  then  a  function  of  /i  and  k.  Neutral  solutions 
then  correspond  to  ftr  —  0  and  the  corresponding  value's  of  h.  B  are  the  neutral  \alues 
of  the  neutral  wavenumber  and  angular  displacement.  In  fart  (^.27)  is  (piite  similar  to 
the  eigenvalue  problem  .solved  by  Seminara  and  Hall  (!97()).  The  latter  authors  were 
conrerned  with  the  stability  of  t*(e  flow  around  a  torsionally  oscillating  cylinder.  The 
eigenvalue  value  problem  /r  =  fi{k\  B)  associated  with  (3.26)  is  itlentiral  to  that  governing 
the  stability  of  a  vertically  oscillating  Roussineset  fluid  of  Prandtl  number  unity  subject 
to  a  time  periodic  temperature  heating  at  the  wall.  If  the  vertical  oscillations  are  replaced 
by  a  steady  gravitational  field  then  we  obtain  the  eigenvalue  problem  discussed  by  Hal! 
(1985).  It  is  of  interest  to  note  that  in  that  case  the  growing  modes  correspond  to 
subharmonic  disturbances. 

A  numerical  investigation  of  the  eigenvalue  problem  (3.27)  showed  that  the  only 
growing  disturbances  have  =  0  .so  that  the  disturbed  flow  is  synchronous  with  tlie 
basic  flow.  Our  calculations  show-ed  that  the  minimum  value  of  B  is  givc'u  by  B  =  2.9 
so  that  at  high  frequencies  the  boundary  between  stability  and  instability  is  given  by 


cv 


2.9 

vr  +  •  •  • 

<I>  4 


\V('  postpone  a  comparison  of  the  low  and  high  frecjuency  predictions  found  above  to  tlu' 
numerical  solutions  of  (3.1)  until  the  next  section. 


4  Numerical  solutions  of  the  eigenvalue  problem 

for  $  =  0(1) 

Oil  the  basis  of  Flotiuet  theory  we  aiiticipale  that  solutions  oi  (.'M  )  may  iie  louiul  in  I  In 
form 

111' 

—  X 

and  the  sign  of  [ir.  the  Floquet  expoiumt.  then  <lelermiiu>s  (he  stability  (iiara<  teiisii(  ol 
the  flow  in  question.  We  obtained  values  of  /i  by  substituting  for  {('A'}  from  (  hi  j  into 

(3.1)  and  solving  the  infinite  si't  of  coupled  ordinary  differential  equatiijus  obtained  by 
equating  like  powers  of  r‘‘  by  a  sliooting  procedure.  Bi'causi'  of  tin'  symnK'tries  of  the 
basic  stale  it  is  possible  to  show  that  the  possible  eigruifunct  ions  {(  nil))-  F)<(,v))  are  eii  her 
odd  or  even  functions  of  y.  This  result  was  used  to  reduce  the  interval  over  which  these 
functions  must  be  calculated  to  [0,^].  However  note  that  all  the  results  we  obtained 
correspond  to  even  modes  in  y  and  the  corres|)onding  Flofiiu't  I'Xpoinuit  was  found  to 
be  purely  real.  The  latter  result  means  that  the  disturbances  are  syuchroiu>us  with  tin' 
basic  state;  we  note  here  that  the  experimental  investigation  of  BM  found  no  evidence  of 
subharmonic  instabilities.  Finally  before  presenting  our  results  \\v  note  that  the  numlx’i 
of  Fourier  terms  used  in  the  truncated  form  of  (4.1)  and  the  number  of  grid  ixiints  in 
the  Runge-Kutta  integration  scheme  were  varied  until  convt'rgence  was  achicwasl. 

In  Figure  (4.1)  we  show  a  sequence  of  neutral  curves  in  the  A'  —  <I>  plane  for  seviual 
values  a.  We  see  that  there  is  a  minimum  value  of  <I>  on  each  neutral  curve.  abo\-e  lhes(> 
curves  exponentially  growung  modes  exist.  If  a  is  varied  we  can  compute  the  o  —  <l>  locus 
of  the  most  dangerous  mode.  This  curve  is  shown  in  Figure  ( 1.2)  and  is  labeled  as  7’!. 
In  this  Figure  we  also  show'  some  of  the  experimental  results  of  B.M. 

The  labeled  /,  II.  Ill,  IV  and  V'"'  were  given  by  BM  and  re])resent  rough  boundaiT's 
between  different  flow  states.  Below  /  no  roll  state  could  be  observed,  whilst  aliove  this 
curve  weak  rolls  could  be  seen  though  their  amplitude  did  not  increase  significantly  until 
II  was  reached.  Regimes  associated  with  weak  rolls  VTA  are  <lenot('d  by  <  ircles  in  f  igure 

(4.2) .  In  fact  we  note  that  the  small  circles  denote  stall's  wlu're  any  horizontal  structure 
was  barely  visible  whilst  the  intermediate  circles  denote  moderate  am|)litude  weak-rolls 
with  defects.  The  large  circles  denote  defect-free  weak  rolls  of  moderate  anqilitmh'.  It 
appears  that  the  theoretical  curve  7’1  prerlicts  the  on.set  of  tlu'  st  rong  roll  cells  assoi  iali-il 
with  II.  Despite  an  exhaustive  search  wecoulo  find  any  amplifying  modes  corresiionding 
to  the  weak  roll  onset  observed  experimentally.  We  note  here  that  no  particular  time 


depeiitU'uce  of  llic  perturhatiou  was  impostnl  in  v)\vv  cal»'vilaliv)iis  su  viia'  ii  siil)iianiiun!i 
or  suporharnioiiic  modes  were  unsf,al>le  they  would  liav<‘  heen  captured  Iw  the  numerical 
scheme.  We  couedude  then  that  thi'  c\irve  /  ol  B.M  is  to  ix’  asstx  iatc'tl  with  (Uid  eifects 
ill  the  experiment,  certainly  the  fact  that  BM  state  that  tht'  roll  am|ditude  d<;es  not 
increase  significantly  until  II  is  cro.ssed  would  tend  to  support  this  conclusion. 

The  diamonds  (o)  in  Figure  (1.2)  denote  the  iiositioii  where  the  niiset  ol  strong 
rolls  was  observed  by  BM  alter  a  slight  inerea.se  in  <I>.  For  o  S  90"  B.M  found  that 
hysteresis  took  place  and  that  the  strong  rolls  did  not  disappear  until  a  lower  \alue 
of  <F  was  achieved.  Tlu'se  points  are  denoted  by  f  in  Figure  (4.2).  The  symbol  II 
was  used  by  13M  to  denote  regions  where  front  propagation  was  observe'd.  here  strong 
rolls  consumed  weak-rolls  as  a  travelling  wave  front  as  4>  was  slowly  increased,  d  he 
symbols  0.5’  were  used  by  BM  to  denote  straight,  defei't  free  rolls  and  straight  rolls 
with  defects  respectively.  In  some  cases  BM  observa-d  wavy  strong  rolls,  tliese  occurred 
in  tlie  same  o  —  region  as  the  straight  rolls  and  these  states  are  denoted  for  dilTerent 
sidewall  conditions  by  A,  IFF'.  At  a  fixed  value  of  o  the  strong  rolls  exhibited  multiple 
superimposed  wavy  modes  [MWV)  which  became  more  disorganized  a,s  <I>  increased. 
Above  curve  IV  strong-rolls  with  superimposed  turbulence  were  observed  by  BM. 

In  Figure  (4.3)  we  compare  the  critical  wavenumbers  of  our  theoretical  predictions 
with  the  observatiotis  of  B.M.  We  see  again  that  the  otiset  of  the  strong  roll  state  again 
correlates  well  with  the  theoretical  work. 

Finally  we  note  that  the  most  dangerous  modes  predicted  by  the  asymptotic  theories 
for  4*  <C  1 ,  4>  1  are  denoted  by  the  curves  A  \  and  .42  respectively  in  Figure  (4.2).  We 

see  that  the  high  frequency  prediction  agrees  well  with  the  finite  calculation  whilst  the 
small  <I>  prediction  is  not  particularly  accurate  at  the  largest  value  of  a  used.  However 
since  the  latter  theory  is  based  on  o  1  and  we  have  computed  only  for  o  S  3  we 
pr('sume  that  tlie  difference  is  because  the  asymptotic,  regime  has  not  yet  been  achieved. 
Indeed  when  cr  ~  3  the  curve  Tl  has  0  ~  .'30  which  means  that  the  unsteady  boundary 
layer  has  a  thickness  of  about  thus  the  quasi-steady  response  of  the  basic  state  is  not 
yet  operational. 

In  Figures  (4.4a,b,c)  vve  show  that  the  first  few’  Fourier  mode's  of  (/,  V  for  the  most 
dangerous  mode  at  a  =  1.1,  1  3,nd  5.  Note  that  these  functions  are  even  about  y  =  0 
and  that  are  zero  when  u  is  an  even  integer. 

In  order  to  see  wliether  the  wavy  and  tiirFudent  states  observa-d  by  BM  an'  related  to 
travelling  wave*  disturbances  we  shall  in  the  next  .section  iliscuss  the  possible  exist('nc<' 


II 


of  Tollmien-Schlichtiiig  instabilities. 


5  Tollmien-Schlichting  wave  disturbances 

Here  we  investigate  the  possibility  that  travelling  wave  (listurl>an(es  are  responsi!)le 
for  the  onset  of  instability  in  the  flow  in  a  flajiping  rectangular  tank.  We  rest  rid  om 
attention  to  two-dimensional  waves  and  therefore  set  d,  —  H'  =  0  in  (2.l'{),  If  the 
pressure  is  eliminated  from  the  x  and  ij  momentum  e(iuation.s  we  find  that  V  satisfies 

W  {aj  _  V  =  {u  +  Wa,}  {ai  _  a^}  v  -  c.i , 

where  we  have  defined  the  Reynolds  number  R  by 

R  —  a4>. 

We  see  that  the  terms  proportional  to  sint  in  (2.13)  do  not  contribute  to  (5.1)  whirli 
is  therefore  the  generalization  of  the  Orr-Sommerfield  eqtiation  to  an  unsteady  parallel 
flow  u  =  u{y,  t) 

The  equation  (5.1)  may  be  solved  using  the  approach  of  Hall  ( lf)78)  who  used  Floquel 
theory  to  convert  the  same  equation  for  a  Stokes  layer  mean  flow  into  an  infinite  sequence 
of  coupled  Orr-Sommerfield  equations.  Here  we  shall  restrict  our  attention  to  solutions 
of  (5.1)  for  small  and  large  values  of  the  frequency  parameter  4>. 

In  the  high  frequency  limit  we  recall  that  u  is  given  by  (2.1)  for  (y  +  j)  =  0(<I>"2).  A 
similar  asymptotic  form  applies  in  the  upper  boundary  layer  so  that  any  instability  will 
be  localized  near  y  =  ±|.  In  fact  u  given  by  (2.10)  is  exactly  the  Stokes  layer  velocity 
profile  for  the  case  when  the  flow  is  driven  by  an  oscillatory  pressure  gradient.  Hence  if 
we  define  rj  =  ^^{y  +  |)  and  let 

A  =  $H,  R  =  <i>+iR 

then  (5.1)  reduces  to 


which  is  to  be  solved  subject  to 

V  =  V,,  =  0,  rj  =  0,oo.  (5.3) 
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rh('  partial  cliirerential  system  (').2)-(r).3)  is  identical  to  that  guveniitiL’,  the  iiistahility  ul 
a  Stokes  layer  to  rollmien-Sclilichting  waves;  see  Kerezek  and  Davis  {  IS71),  Hall  I  197S). 

Thus  the  stability  of  tlie  (low  in  a  rectangular  tank  iiapjiing  at  high  frtn|nenrics  i'- 
governetl  by  the  equations  which  determine  the  stability  of  a  Stokes  layer.  The  Flofjuet 
analy.sis  of  (2.14)  given  by  Hall  (1978)  suggests  that  a  Stokes  layt'r  is  stalile.  on  the 
other  hand  the  quasi-steady  ai)praarhe.s  of  Kerczek  and  Davis  (1971)  show  that  instan¬ 
taneous  profiles  can  !)e  highly  unstable  because  of  their  inlh'-xiiuial  nature,  ’{'he  results 
of  tliese  different  approaches  can  be  reconciled  by  noting  that  the  <iuasi-steady  siilutions 
cannot  lie  connected  to  the  Floquet  solutions  l)y  extending  them  over  a  whole  period. 
Nevertheless  the  results  of  the  quasi-steady  calculations  are  consistent  with  experimen¬ 
tal  observations  and  suggest  instability  will  occur  for  part  of  the  [jeriod  whenever  H  is 
greater  than  about  200.  This  suggests  that  at  high  fretpiencies  localized  inst  altilities  iis 
the  form  of  Tollmien  .Schlichting  wavc^  will  occur  when 

200 


a  > 


whilst  roll  modes  occur  when 


2.9 

t>  >  — r- 
4)4 


It  follows  that,  in  an  experiment  with  4>  fixed,  transition  will  probably  be  cau-seHl  by 
Tollmien-Schlichting  waves  for  small  enough  values  of  cv.  However  o  must  be  less  than 
about  .05  for  the  Tollmien-Schlichting  wave  to  become  dominant;  this  regime  was  not 
investigated  experimentally  so  it  is  not  surprising  that  the  asymtotic  prediction  given 
above  is  off  the  scale  of  Figure  (4.1). 

Now  let  us  turn  to  the  low  frequency  limit  4>  — )  0,  in  this  limit  u  is  given  by  (2.1 1) 
and  (5.1)  may  be  written  in  the  form 


1 


iXR 


~  D  —  cost  ^ j-  F  —  costllyy  I'  — 


where  u  =  1  —  4?/^}. 


/?  =  4>/7. 


0.  (5.4) 


(5.5) 


The  equation  (5.4)  is  to  be  solved  subject  to  the  conditions  that  17  should  vanish 
at  y  —  ±^.  The  slow  time  variation  of  the  basic  state  can  be  taken  care  of  for  the 
disturbance  by  a  W KB  approach,  we  therefore  expand  V’  in  the  form 

V  -  [Vo(y,  0  -p  4)Vj(y.  0  -f  .  .  . 

and  r(t)  is  then  fletermined  by  the  instantaneous  Orr-Sommerfic'id  eigenvalue  ))robl('m 


?XB 


{yj  A-q'  c,  =  ),■<»(  -  0  {g'  -  A'd 


rosiUyyWt  -  0, 


Hi 


V„  =  v;,^  =  u.  //  =  ().,  1. 

riius  llu' effective  Revtiolcis  imniix'r  associatoil  with  tli<' xf'hait  y  iiehl  (7{//)  is  In 

Figure  (5.1)  \V(' siiovv  the  iKuitral  curve  A  =  A{/fces/)  which  tiiarks  th(“  luiunihti  x  hciwrcn 
locally  growing  aiul  decaying  solutions.  I'hus  whenever  t  and  U  are  su<  h  that  l/Vcos/:  i> 
greater  than  its  valiK'  on  the  neutral  curve  at  a  lixe<|  A  the  solution  is  ha  ally  growing. 
We  further  note  that  Figure  (5.1)  is  the  neutral  curvo  fcjr  the  cur\'e  for  tin-  lolhnien 
Schlichting  wave  instability  of  the  flow  a  chaniud  r»)tating  with  a  steady  angular  veloi  ity. 
In  addition  it  should  be  notetl  that  all  the  eigenvalues  we  found  had  zercj  wa\ espc'cil  >o 
that  the  instability  wave  corresponds  to  a  standing  wave  instability,  liowesiu-  huge  we 
choose  R  it  will  be  the  case  that  cos/|  is  sufficiently  small  for  tin*  part  v>f  flu-  cycle  so 
that  the  disturbance  is  locally  decaying.  The  integrated  value  of  lh('  growth  rate  over 
a  cycle  then  determines  the  stability  property  of  tlu'  flow  according  to  I'lo<(uet  thetu>-. 
Note  that  it  is  sufficient  for  us  to  consider  only  ^  a  ])eriod  of  the  basic  flow  so  that  the 
neutral  solutions  based  on  Flociuet  theory  are  given  by 


(5.7) 


where  /{}  denotes  the  imaginary  part  of  a  complex  (juantity.  In  Figure  (5. 1 )  we  also  show 
the  neutral  curve  obtained  by  the  imposition  of  this  condition.  W’c'  see  that  instabilil\' 
occurs  for  R  >  '1965.  and  unstable  morles  occur  over  a  fiidte  range  of  value's  i)f  llu' 
wavenumber  A.  The  fact  that  the  band  of  unstable  wavt'iiumbers  is  finite  is  a  din'd 
consequence  of  tfie  fact  that  17  has  ati  inflection  point  .so  that  at  any  instant  in  lime 
when  cost  7^  0  at  sufficiently  high  values  of  R  the  instantaneous  lu'utral  problem  has  a 
mode  with  Cr  =  0.  Finally  we  note  that  our  calculations  predict  tlu  onset  of  I'ollmie'm 
Schlichting  instabilities  when 


cv  > 


3965 


(5.S1 


The  unstable  region  predicted  by  (5.8)  is  off  the  scale  in  Figure  (1.2)  so  it  would  api>ear 
that  in  BM  the  wavy  and  turbulent  regimes  are  not  associated  with  rollmif'ii-Schlif  lit ing 
waves.  However  it  should  be  remembered  that  (5.8)  is  valid  onlv  for  <t>  1 ,  at  0(1) 

values  of  the  frequency  parameter  the  stability  of  tlx'  basic  state  can  (uily  Ix'  deter¬ 
mined  Iry  solving  (5.1)  numerically.  VVe  do  not  pursue  that  calculation  here'  since  our 
a.symptotic  results  suggest  that  Tolhiuen-Schlichtiiig  waves  ar<'  not  important  in  the  ex 
perimental  ranges  investigated  Iry  HM.  Nevertheless  the  FfiKpiet  solutions  of  (5.1  )  ar<'  of 
some  interest  since  we  know  that  at  small  values  of  unstable  solutions  exist  whilst  at 
large  <I>  the  ef|uation  (5.1)  g(;V(’rns  tin*  instability  of  a.  Stok('s  lay<'r.  Since  no  unstable 
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Floqurt  solutions  have  ever  heeii  obtained  for  the*  laMt'r  piuljh'in  it  '.vill  be  of  inleiot  to 
determine  here  tlie  small  <I>  unstalile  solutions  ixsome  stabilized  at  larger  <h. 

6  Conclusions 

Our  investigation  has  shown  that  the  flowing  an  oseillating  fluid  tank  is  sus(  ('[)t  ibie  to 
at  least  two  types  of  instability.  The  first  moile  is  the  roll  mode*  haxing  (fll  bcamdaries 
parallel  to  th  x  —  y  plane  whilst  tlie  other  instability,  the  wave  morle,  is  pericxlic  in  the  s 
direction.  Furthermore  in  the  low  fretjiiency  limit  the  wave  mode  is  stationary  so  that  t  he 
instability  takes  the  form  of  rolls  which  are  now  parallel  to  the  y  —  :  plain'.  I'lie  onset  of 
the  strong  mode  observed  by  BM  is  explained  by  the  most  ilangerous  linear  disturbance 
discussed  in  Section  4  of  this  paper.  We  believe  that  the  weak  roll  observc'd  by  H.M 
is  a  manifestation  of  end-effects  in  their  apparatus  and  is  thert'forc'  not  accessible  to  a 
linear  instability  analysis.  The  roll  modes  vve  have  discus.sed  have  a  close  relat  ionshij) 
with  centrifugal  and  cotivective  instabilities  in  titne  periodic  boundary  layers  and  it  is  of 
interest  to  determine  the  destabilizing  mechanism  in  the  present  situation.  In  fact  we  see 
iti  (3.1 )  that  the  terms  on  the  right  hand  side  of  the  F  equation,  whici,  are  responsible  for 
the  instability,  arise  from  the  Coriolis  terms  in  the  N'avier  Stokes  equation  written  down 
in  the  rotating  frame.  Thus  the  roll  mode  is  produced  by  Coriolis  effects.  On  the  other 
hand  the  wave  mode  is  associated  with  an  inflection  point  instability  in  both  the  high 
and  low  frequency  limits.  More  precisely  the  wave  instability  discus.sed  in  §•')  at  small 
values  of  is  an  inviscid  instability  associated  with  the  inflectional  velocity  |)rofiIe  U{y). 
Our  calculations  suggest  that,  low  frequencies,  this  mode  does  not  play  a  significant 
role  in  the  highly  nonlinear  stages  investigated  experimentally  by  BM.  However  the 
wave  instability  might  be  more  unstable  at  finite  values  of  the  frequency  paramet<'r.  this 
possibility  has  not  been  investigated  numerically  but  certainly  the  results  of  B.M  suggest 
this  as  a  strong  possibility.  At  high  frequencies  the  wave  mode  is  a  locally  unstable 
Stokes  layer  instability,  BM  do  not  give  any  results  which  suggest  that  this  instability 
is  present  in  their  experiments.  We  note  that  this  is  a  highly  localized  mode  and  so  if  it 
were  present  in  the  experiments  it  would  almost  certainly  be  detected. 

In  addition  to  the  higher  order  linear  roll  modes  and  travelling  wave  disturbances 
there  are  other  candidates  for  the  .secomiary  and  tertiary  instabilities  obsc'rved  in  the 
exjjeriments.  We  refer  to  the  inviscid  modes  induced  liy  finite  am|>litude  vortices  in 
boundary  layer  flows.  The.se  modes,  which  were  <liscussed  by  Hal!  and  Horseman  (1901), 
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arise  when  a  amplitude  vori<'x  modifies  tiie  underiying  houndarv  la\<'r  so  as 

make  it,  imstahle  to  Rayleigh  waves,  d  iu*  s<*eoiidary  instahility  in  that  proldem  tfnd>  to 
he  treated  at  particular  spanwist'  locations  hut  tliere  is  nt)  suggestion  that  this  f,vi>e  of 
localization  occurs  in  the  fla])ping  tank  prtddem. 
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